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On the Diffraction and Reflection of Waves and 
Pulses by Wedges and Corners 

F. Oberhettinger l 

Various problems arising in the theory of the excitation of a perfectly reflecting wedge 
or corner by a plane, cylindrical, or spherical wave, are dealt with. The incident wave is 
represented by a line source (acoustic or electromagnetic) parallel to the edge. The spherical 
wave is emitted by an acoustic point source or by a Hertz dipole with its axis parallel to the 
edge. The case of an incident plane wave field is obtained as the limiting case (for large dis- 
tances of the source from the edge) of the cylindrical or spherical wave excitation. 

1. Introduction 

Sommerf eld's classical problem of the excitation of a perfectly reflecting half plane by a 
plane (acoustic or electromagnetic) wave field is based upon the theory of multivalued solutions 
of the wave equation [34]. 2 Subsequently this method was extended to the problem of the dif- 
fraction of cylindrical or spherical waves by a half plane [3] or by a wedge [44] (in the latter 
case also for a plane incident wave). Here the incident cylindrical wave is emitted by a line 
source parallel to the edge (of the half plane 4 or wedge) and the spherical wave is emitted by a 
point source (in the acoustic case) or by a Hertz dipole (in the electromagnetic case) parallel 
to the edge. The solutions for the half plane problem are given in the form of integral expres- 
sions involving a real variable, while the corresponding wedge problem solutions are given as 
contour integrals. 

The problems described above have been treated ever since with the aid of various techniques. 
The wedge problem as treated by MacDonald [22] requires the conventional methods of nine- 
teenth century mathematical physics only. Here the solution is given by an infinite series and 
this is transformed into Wiegrefe's complex integral. For the special case of a half plane (angle 
of the wedge equal to 2 tt), reductions to integrals involving a real variable are given. Another 
approach is a generalization of the method of the images [1, 16, 43]. 

At first the solution for a wedge of an angle w/m with positive integer m is established as 
the sum of the effects of the original source and its 2 m-\ images. This sum is replaced by a 
contour integral which is extended so as to obtain a formula valid for any positive m and hence- 
forth for any angle a of the wedge. Again, the solution is given by a contour integral. In 
more recent times a number of attempts have been made to apply integral transform methods 
[20, 17] and integral equation methods [23, 7, 14, 13] to the half plane diffraction problem. 
Finally, the problem of the diffraction by a perfectly conducting half plane of electromagnetic 
waves emitted by an arbitrarily oriented electric or magnetic dipole has also been solved [32, 
39, 15]. 

Sommerfeld's solution (plane waves incident on a half plane) is given in a form readily 
available for numerical computation inasmuch as it contains only the tabulated error func- 
tion besides elementary functions. The solution of the other cases can either- be given in the 
form of infinite series (usually suitable for the numerical evaluation of the near field) or in the 
form of integral expressions. A reduction of the latter to an asymptotic series has to be carried 
out if approximate expressions are desired for distances that are large compared to the wave 
length. Expressions of this form have been derived for the case of the diffraction of a plane 
wave by a wedge [44, 29]. (An asymptotic expansion which fails to give the proper behavior 
in the neighborhood of the boundaries of geometric optics is ignored here [36].) The main 
feature of these asymptotic series is that their terms are generalized error functions. The work 
referred to above is concerned only with time harmonic problems. 
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A second group of investigations centers around problems in which the dependence upon 
time is not periodic. This is the case of the diffraction of a pulse or a transient wave by a wedge 
or a half plane. This group of investigations was initiated by Sommerfeld [35]. The following 
contributions are listed here: Normal incidence of a plane pulse on a half plane [21,11]; two- 
dimensional pulse incident on a half plane [12]; plane pulse (unit step function) incident on a 
wedge or a corner [19,24, 25]; arbitrary pulse incident on a wedge or a corner [18]; cylindrical 
Dirac pulse incident on a half plane [38]; and spherical pulse (unit step function) incident on a 
half plane [40]. This paper will deal with various problems arising in the theory of the exci- 
tation of a perfectly reflecting wedge or corner by a plane, cylindrical, or spherical wave field. 

The incident cylindrical wave is represented by an acoustic or electromagnetic line source 
parallel to the edge. The spherical wave is emitted by an acoustic point source or by a Hertz 
dipole with its axis parallel to the edge. The case of an incident plane wave field is obtained as 
the limiting case (for large distances of the source from the edge) of the cylindrical or spherical 
wave excitation. 

In section 2 of this paper, various representations for the time harmonic solution of the 
problems outlined before are given with special emphasis to a representation upon which the 
investigations of sections 3 to 6 are based. A straightforward method is applied based upon the 
representation of the incident field in the form of a Kantorovich-Lebedev transform. 

Section 3 contains a detailed discussion of the results obtained in section 2. Relations to 
geometric optics are discussed. 

Section 4 contains an asymptotic expansion for the far field excited by the incidence of a 
plane wave on a wedge. This expansion is not of the Wiegrefe-Pauli type inasmuch as it con- 
sists of two error function terms and an asymptotic series with inverse powers of the variable. 
Remarkable is the fact, that this series represents asymptotically (for large distances from the 
edge) the deviation of the wedge field from the corresponding half -plane field. 

Section 5 contains expressions for the energy radiated from a Hertz dipole in the presence 
of a perfectly reflecting wedge or corner (antenna and corner reflector). 

Section 6 gives the pulse or transient solutions corresponding to the time harmonic prob- 
lems investigated in section 2. It is shown that previously known results can be obtained as 
special cases of the expressions given here. 

2. Time Harmonic Wedge Diffraction Problems 

The case of the excitation of a perfectly reflecting wedge by a (acoustic or electromagnetic) 
line source parallel to the edge of the wedge is considered first. (The solution for the excitation 
by a point source can then be obtained by a further integration.) C.ylindrical coordinates 
(p,<t>,z) are introduced such that the z axis coincides with the edge of the wedge and the walls of 
the wedge are bounded by <£=0 and 4>=a. Let the locations of the line source and the point 
of observation be characterized by Q(p',4>') and P(p,<l>), respectively (fig. 1). Let the field 
inside the wedge be periodic'of^the form 

<S>=U(p,4>)exp(ia>t), (1) 



<p=a 




*P=0 



Figure 1. 
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where U satisfies the two-dimensional wave equation 

52/7 1 XT J 1 52f7 

AU+k > U== *U + l ^ + J_ |^ W=0 . (2) 

A solution of (2) representing divergent waves is of the form 

U=ex V (±iv<i>)FLV>(kp). (3) 

H v {2) (kp) is the second Hankel function, and v is an arbitrary parameter. Following a pro- 
cedure used elsewhere [26], the substitution 

k = 27r/\=w/c=—iy (4) 

is introduced. This means that instead of (2) and (1) 

AU-y 2 U=Q, $= L r (p,0)exp( 7 cO (5) 

is considered with a particular solution corresponding to (3) 

U=exp(±iv<f>)K v (yp), (6) 

where 

K v (z) = -±iTroxv(-i\v<K)H v u\-iz) (7) 

denotes the modified Hankel function. The transition from k to 7 amounts to the transition 
from a wave problem to an exponential decay problem. At first the exponential decay prob- 
lem is solved, and the wave problem solution is established by returning from 7 to k. This 
will not only simplify questions with respect to convergence but the investigations in sections 
4 and 6 are based upon the exponential decay solution and not on the time harmonic solution. 
The incident field in P(p,<t>) due to the line source <?(p',0') is given by 

[/'=#o (2) (fcr)=-#oMp 2 +P ,2 -2pp' cos {4>-4.')Y*}- (8) 

7T 

Let the total field U be represented as a sum of the incident field W and a reflected 
field U T , where U r has no singularities within the wedge: 

U=U t +U r . (9) 

In order to represent U l in a form suitable for the problem considered here it is advisable 
to express the modified Hankel function on the r. h. s. of (8) not by means of its well-known 
addition theorem [9, p. 44] but by the equivalent expression ((3A) appendix): 

U'=- Zoh[p 2 +p' 2 -2pp' cos (4>-</>')] 1/2 } 

TV 

= Ut-* £ K iv (yp)K iv (7p') cosh [v{ir-\<j>-cj> f \)\dv y (10) 

which is of the form (6). Let the reflected field be represented in the same form: 

U'=4iir-* ^^.(T^^^TpOI/iWexpM+^WexpC--^)]^. (11) 

The so far unknown functions j Y and/ 2 are to be determined by the boundary conditions. 

There are two cases of boundary conditions: (a) U=0 and (b) dU/d<l>=0 on the walk 

</>— O^and <j>=a of the wedge. With electromagnetic waves on a perfect conducting wedge, 
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case (a) refers to the electric field, case (b) refers to the magnetic field. With acoustic waves, 
case (a) refers to free wall, and case (b) refers to rigid wall. 
This leads to 

fi+f2^-coah[v(ir-4>')] 

J 1 exp(va)^rf 2 exp( — vot) = — cosh[^(x— a + 0')] (12) 

and 

yi-y^-sinhWT-*')] 

f t exp(va) — y 2 exp '(—va)= sinh[^(ir— a + <£')]. (13) 

for cases (a) and (b), respectively. It may be noted that the boundary conditions have been 
satisfied by putting the integrand equal to zero. This can be justified by verifying that the 
solution thus obtained satisfies all necessary conditions [31]. But it may also be remarked 
that ?7is represented as a Kontorovich-Lebedev transform [9, p. 75]. After j x and/ 2 have been 
determined by (12) and (13) the solution is then, by (9), (10), and (11), 

U=±iir- 2 C K^p)Ki V (yp' f ml x ^ {coshfKa-l^-^DlTcoshfKa-^-^)]^^- (14) 

Jo si nil [av) 

The upper and lower sign is valid for cases (a) and (b) 7 respectively. The relation ((9A) 
appendix) shows that the convergence of the integral (14) at its upper limit for a real and 
positive 7 is secured when 

and this region contains the whole region of the wedge. If one returns now to the wave problem 
replacing y by ik according to (4) one arrives at the result 

U=-i Pexp {Trv)H^{kp)m v ) {k P F ) S ! n l 1 ( /^ {cosh Ka-^-^'DlTcosh [v(a-<l>-<f)')]}dv. (15; 
Jo sinn {av) 

The behavior of the integrand at its upper limit leads by [9, p. 88, form. 22] to the condition 

x^ 10-F <t>'\ = 2a-|-7r. 

This restricts the validity of (15) to the case of a wedge of an angle a larger than w with 
the point of observation P in the geometric shadow of the source Q. The expression (14) 
shall now be transformed into an infinite series. For this purpose the elementary relation 
[2, p. 368] 

coshKa-fl)] v " € n cos(nire/a) , 

sinh (av) at^s n 2 + (mrja) 2 

valid for 0^0^ 2a is used. Term by term integration yields 

£K t ,(y P )K tv (y P ') ™^"~f ) ] sinh (Tv)dv=\ ±e n cos (nrd/a) £ p -^^K ir (yp')K ir (y P )dv. 

The integral is known f 10, vol. 2, p. 176]. (Note that this formula is not only valid for positive 
integers but for Re n^>0). Hence 

I **+(! i V l K iv (yp)K^yp , )(lv=lTr 2 I nlr/a (yp)K n7r/a (yp f ) 

Jo v'-tCnir/a)' Z 

for p<Cp' and for p> p the same formula with p and p interchanged. Therefore, the field U can 
be written as 

^=f J^f 008 ^^ - ^* 008 ^ (*+*)']} I».,*(yp)K„„(yp') (16) 



for p< p' and for p > p the same formula with p and p' interchanged. The expression (16) is 
valid for arbitrarily complex 7 and if 7 is replaced by ik the solution for the corresponding wave 
problem becomes 

k T =^g«»/oos[^^ (17) 

for p<Cp' < and for p > p r the same expression is used with p and p interchanged. The preceding 
expressions (16) and (17) represent rapidly convergent series if the arguments of the Bessel 
functions are small. The expressions (16) and (17) shall now be transferred back into another 
integral expression which is the basis of the investigations in sections 3 and 6. For this purpose 
the product I„ir/a(yp)K nr / a (yp') in (16) is replaced by an integral representation ((1A) appendix). 
This leads to the following expressions. 

CO 

«=S 6 « cos (nir8/a)I nv/a (7p)K nir/a (yp') 

11 = 

= *" 1 2«ii cos M/o) f 7l [7(p 2 + p' 2 -2pp' cos ./•)' 2 ] cos (mrx/a)dx 

// = () Jo 

-^S «» COS (/<ir0/a) sin (ri7r 2 /a) | exp (— nrx/a)K Q [y(p*+p'*+2pp' cosh .r)^ 2 ]r/.r. (18) 
>/=() Jo 

The second sum vanishes when a=irjm i m=l, 2, 3, .... The two sinus on the r. h. s. 
of (18) are denoted by .*, and s 2 . 

Interchanging the order of summation and integration one finds 

s l =(2v)' 1 lim j 7v [7(p 2 +p' 2 -2pp' cos r) 1 / 2 ] £> n /cosfti | (a +0)1+ cos f/,^ (.,-0) 1 \ dx, 



s-o^-tt" J CC /v [t(p 2 +p /2 +2pp / cosh ./-) I/2 1 

• X)exp(— /i7r.//a)^ sin n- (tt+0) +sin n- (t—6) >dx. 

The sum in the integrand of (20) can be summed [2, p. 203] and the result is 

s 2 =-(2t)-i J V [y(p 2 +p' 2 +2pp' cosh x) 1 *] 

sin [~^ (a— 0) ~] sin [~^ (tt+ 0) 1 

cosh (tx/o) — cos (71—- 0) cosh (irx/a) — cos - (7r+0) 



(19) 



(20) 



-<*r. 



(21) 



s, as given by (19) can be expressed as a sum of a finite number of terms. The result is (appen- 
dix) 



, 1 = a7r- i J^KJ'r[p i +p''-2pp'cos(2ar+ej}"A- 
The limits of the summation in (22) are 

r, 1 



(22) 



where the well-known symbol 9 means the largest positive or negative number which 
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is smaller or equal to —= — The sum (22) is zero when r{^>r 2 and the first or last term of the 

sum has to be halved when r x or r 2 is an integer. With (18) to (23) the expression (16) for the 
field can also be written as 

U=F(4>-4>')TF(4>+<I>'), (24) 

with 

F{8) = 2iTr-' i^^ { 7 [p 2 +p /2 - P p / cos(2ar+^l 1 / 2 }-^- 1 V ' K {) [y(p 2 +p' 2 +2pp' cosh x) 1 '*] 

r = n JO 



sinf^-fl)! 



sm 



[i (,+»,] 



cosh (irJc/a)—cos - {ir—d) eosh (wx/a) —cos - (ir + 6) 



dx, 



(25) 



and ri and r 2 defined by (23). 

2.1. Spherical and Plane Field Excitation 
The formula [9, p. 95, formula 53] 

( r 2 +3 2)-i/2 exp [_ T ( r 2 +2 2)i/2]_27r- 1 K [r (y 2 + v 2 ) 1/2 ] cos {zv)dv, 



(26) 



expresses the generation of a spherical wave function by means of a cylindrical wave function. 
To obtain the field U corresponding to an incident spherical field multiply the r. h. s. of (8) 
and (25) by cos (zv), replace y by (y 2 +# 2 ) 1/2 and integrate the resulting expression with respect to 
v from to oo . This operation performed on (8) gives, using (26), 



f* <x> 

217T- 1 K [r(y*+v 2 yi i ]cos(zv)dv=i (r 2 +2 2 )-i e .xp[-Y(r ! +2 2 ) 1/2 ] 



(27) 



i. e., an incident spherical field originating from a source point Q(p' , </>' ', 0). The same operation 
performed on (25) gives similar expressions. (For similar procedures, see [42, 5]). It follows 
then that in the case of an incident plane, cylindrical or spherical field the field U within the 
wedge of an angle a can be written in the form 



rwX0-<w^W>+</>') 

F(6) is expressed as the sum of two functions F x and F 2 : 

F(6) = F 1 (d) + F 2 (d), 



(28) 



where 



F 2 (e)=-(2a)- l j o ° > H(x) 



(29) 

Fi(6)=f:f(e+2ctr), (30) 

in f^ (tt-0)~] sinf- (t+6)~] 



sm 



cosh (tx/o)— cos - (tt—6) cosh.(irx/a)— cos - (x-ffl) 



»dx. (31) 



The representations (28) are also valid for complex values of y with Re y^O and conse- 
quently for y=ik with ImHO. The functions/ and H in (30) and (31) depend on the char- 
acter of the excitation and are listed below. (The incident plane field case is obtained as the 
limiting case of the incident spherical field case for p'->oo .) 
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Functions j and H for the exponential decay probli m 



(incident field) 


//(,-, 


expfrpcoste— <j>')] 
(plane field) 


exp( — 7pcosh x) 


2iir- 1 K (yr) 

(cylindrical field) 


2/tt- 1 K [y (p 2 + p" 1 + 2pp'cosh X)" 2 ] 


( r 2 _|_ 2 2)-l/2 exp [ _ y ( r 2 + ,2) I ft] 

(spherical field) 


(p 2 + p ,2 + ^4-2 P p'cosli as)-i/i 
•exp[— 7(p 2 +p /2 + 2 2 -f 2pp'cosh .r) 1 -]. 



Functions f and H for the wave problem 



/Op-*') 

(incidenl field) 


// ,- 


e\|)[/A'pCOS(<> -</>')] 
(plane wave) 


exp( /A'pcosh a?) 


Hl> 2) (kr) 

(cylindrical wave) 


// ( ; 2) [A-(p 2 + p' 2 + 2pp'c()sh.r)'-'l 


( r 2 + 2j2)-i/2 ex p[-zfc(r 2 + z 2 ) 1 / 2 ] 
(spherical wave) 


(p 2 + p' 2 +* 2 + 2pp'cosh a?)" 1 / 2 
•exp[ *'fc(p 2 + p' 2 + 2pp'cosh .r • :-)'-!. 



r=[p*+p' 2 -2pp'<-os (*-*')]' 



(32) 



The function /(</>— 0' ) represents the analytic character of the incident plane, cylindrical, or 
spherical field. The upper or lower sign in (28) has to be taken according as an incident elec- 
tromagnetic or acoustic field is considered. The total field U is now represented in three dif- 
ferent forms, (14), (16), (28) and (15), (17), (28) for the exponential decay and the wave prob- 
lem, respectively. Since in the following, (28) will be dealt with exclusively, the expressions 
for U of the form (14), (16), (15), and (17) corresponding to an incident plane or spherical 
field are not listed here. They are obtained in the same manner. The field as represented by 
(28) consists of the sum of two functions F x and F 2 . Each term r=0 of F } represents a field of 
the same character as the incident field (r=0) but originating from one of the images of the 
source Q with respect to the walls of the wedge. The number of the images and their location 
depend on r x and r 2 and henceforth on the configuration. F 2 is given by an integral expression 
which vanishes when the angle of the wedge is of the form ir/m, m=l, 2, 3, . . . . It is easily 
seen that in this case ( r =7<\ (</>—(£' )±Fi(</> T</>') is composed of the effect of the original source 
and its 2 m — 1 images. In case of an incident plane field and an angle 2x (half plane), F 2 can 
be reduced to the error integral (39). 

3. Discussion of F ] and F, for Particular Cases 



It is obvious that the number of terms of the sum (30) can be made arbitrarily large if 
a is arbitrarily small. The case 0<a<V shall from now on be referred to as a corner. It will 
be shown now that for a wedge (7r<a<27r) the sum (30) consists of at most one term. Two 
cases must be distinguished here: case a, a>7r, 0<<p'<C a— 7r » figure 2; and case b, aT>7r, 
a— 7r<<p'<^a:, figure 3. The restriction <t><i\a means no loss of generality. 
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<P=Q-7T 



<p~~tT-y> 




^--7T + <p 



FlCxURE 2. 



<p = 



<^ = a 



ip-.a-7r 



y>=7r-y> 




V>--0 



Figure 3. 



These cases represent the two possible configurations. 

Care a. The direction </>' of the source Q is such, that a geometric shadow region exists. 
The boundary of the geometric reflection is <t> 1 = w— <£' and the boundary of the geometric 
si lado w is <£ 2 = tt + </> ' . 

Case b. The direction 0' of the source # is such that no geometric shadow region exists. 
The boundaries of the geometric reflection with respect to the walls </> = and </> = a are (f)i = 7r—(f) f 
and <l)2 = 2a—T—a / , respectively. 

Let the region rg</>rga be subdivided into three regions I, II, III bounded by the walls and 
the boundaries of geometric optics <£ = $i and <j> = <t> 2 . They are defined by: 

Region I. Og^^ 

Region II. <j> x ^<t>^<j> 2 
Region III. 02^<£^a 



and 



4> l = 7T—(f) / , </>o:=7r + <£' case a, 
cj) l = ir — <£', <t> 2 = 2a—TT—<f)' case 1). 
In these cases the limits r x , r 2 of (30) are easily found to be as listed below. 





Case 


a 




Case b 


Region 


I 


II 


III 


I 


II 


III 


r u r 2 


0,0 


0,0 


0,-1 


0,0 


0,0 


0,0 


run 


0,0 


0,-1 


0,-1 


0,0 


0,-1 


-1,-1 



(33) 

(34) 



The upper or lower line for (/*i,r 2 ) is valid according as d=4>—4>' or d—4>-\-4>' . With these 
results one obtains by (30) for both eases a and b, 



F 1 (4>-<t>')^l'\(<t>+4>')=.f(4>-<t>')Tf(4> + <t>') valid in I, 
FM-<t>')TFM J r4>') = f(<t>-<l>') valid in II, 



(35) 
(36) 
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while in region 1 1 1 



^,(*-0 / ) : F/ T i(* + * / ) = O case a 

F l (<l>-<t>')TF l (<l>+<t>')=f(<j>-<t>')Tf(<l> + <i>'-2a) case b. 



(37) 
(38) 



Equations (35) to (38) together with (31) give the representation for the field Ubj (28) 
and (29) in agreement with the results obtained elsewhere [26]. The ease of a source within 
a corner (0<^a<^w) (fig. 4) and referred to as case c is treated the same way. At first i\ and 
r 2 are determined by (23) and then F x by (30). Together with F 2 by (31 ), the field J r is known 
by (28). 



(£> = a 



V>=<A 




y>=v>o 



*</> = 



Figure I. 



Analogous to the cases a and l> there exist two directions </>, and </>•_> in which a plane wave 
field incident on the corner emerges after multiple reflections on the walls. They are given in 
(50) and (51). 

The second function F 2 (0) given by the integral (31) can be given explicitly when « = 27r 
(half plain 4 ) and when the incident field is a plane field. Since this result will be needed Later 
it is represented again although known [34]. Now, for a plane incident field, /(</>—</>') 
= exp [yp cos ((/>— </>')] and H(x) = exp(—yp cosh x). For <x=2t, 

F 2 (6)= -7T" 1 cos (§0) J cosh ( - xj (cosh ./+cos d)~ ] exp ( — yp cosh x)dx, 

by (31). The substitution cosh x=l-\-t gives 

F 2 (6) = -tt- i 2- 1/2 cos (~$\ exp (-yp) ( r 1/2 (l + cos + O -1 exp (-ypt)dt. 

This integral is known [10, vol. 1, p. 136]. Thus 

F 2 (S)= -| cos (J 0) | sec (J 0)|Erfc [(2 7P ) 1/2 | cos | 0|] exp ( 7P cos 0), 



where 



Erfc (z) = 2ir~ l/2 [ exp (-r 2 )^. 



(39) 
(40) 



Since the half plane case is represented by case a, as described in figure 1, the matching 
expression for /'', in the three regions is known by (35) to (37), and the field U (28) excited by 
an incident plane field on a perfectly reflecting half plane can be written down explicitly. It 
can easily be verified that U can be represented by an expression which is valid over the whole 
range of the half plane observing that by (40) 
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Erf (-z)=2-Erfc (s), 
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(41) 



and so U becomes 

t7=exp [ yp cos fo-0')]{i-J Erfc [(2 7P ) 1/2 cos (fa-**')]} 

Texp [ 7P cos (0+^)]{l-i Erfc [(2 7 p) 1/2 cos(± 0+j4>')]}, (42) 
for the exponential decay problem. The solution for the wave problem is obtained as 
C7=exp [ikp cos (^>-^)]{l-i Erfc [(2ikp) 1/2 cos (§ <$>-\ </>')]} 

=Fexp [ikp cos (<£+</>')]{ 1-i Erfc [(2ikp)^ cos (J *+**')]}• (43) 
(For tables of the error function for complex argument see [6].) 

In general case, F 2 consists of a sum of two integrals of the form 

F 2 (e) = I(a,T-0)+I(a,T+6), (44) 

with 

/(«,#) = -(2a)- 1 I H(x) sin (7^/0) [(cosh (<jrx/a)-eos (w^/od^dx. (45) 

Obviously the denominator of the integrand (45) is zero at its lower limit x=0 when 

*=^«=»±(*±* / )=2wa J w = 0,±l,±2, . . . . (46) 

But the numerator of the integrand vanishes likewise for the same value. It is easily shown 
[27] that 

lim/(a,0=TiH(O) (47) 

and the behavior of (45) is the same for \f/— >±2ria. It may be noted that (47) does not depend 
on the angle of the wedge. For any configuration is always 

0<«+* / <2a, 0<l*-*'|<a. (48) 

This condition restricts (46) to onty two possible directions <j> x and <j> 2 which define the 
boundaries of geometric optics as outlined before. They are defined by 

case a: ^=^—^',(71=0); <fe='H-<£',('rt = 0), 

caseb: <i> 1 = t— <£', (n=0); </> 2 =2a— 7r— <£', (n = l), 

in agreement with (36). A choice of values for the parameter n other than given above is 
incompatible with the configurations under consideration. In order to treat case c it is 
practical to subdivide this case according as the angle a of the corner can be described as 



(49) 



a=T/(2l + e), 1=1, 2, 3, . . . case c a 
a=w/(2l+l + e), 1=0, 1, 2, . . . case c 2 , 
where € is a number between and 1. Then the relation (46) leads to 

<h=\ea-<l>'\ 1 

Vcase Ci (50) 

2 =Min [a (2 -e)-</>' ,«' + €«] 

and 

&=|«(i-*WI 

case c 2 (51) 

<t> 2 = \lm[a(l+e)-<t>',a(\-e) + <l>'] 
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and e is connected with a by (49). In case e=0, this lends to pi = <p2 and particularly, 

<t>i=<f>2=<t> f when a=w/(2l) 

( f )l = (f >2 =^a—(l) / when a=x/(2Z+l). 

An elementary expression for the upper limit of I(a$) can he given. Integration of (45) by 
parts gives 

00 

—TrI{a^/)=H{.r) arctan [cotan (J x^/a) tanh (J x.r/a)] 

o 

— I //'Or) arc tan [cotan (4 x^/a) tanh (ix.r/a)]rfz. 
Since, with the exception of the plane wave case, H(x)->0 when x-^ m it becomes 
I{a\\p) =7T _1 I H' (./) arctan cotan ( ^mf/fa ) tanh ( - wx/a f dx, 

! I(cn4f) | ^ tt" 1 j \H' (.r) 1 1 arctan fcotan Q w^/a) tanh (i tx/o\ Ida :, 



(52) 



and hence 



and, after replacing tanh (^Tr^/a) by its upper limit 1, 

|7(a ] ^)|^|//(0)||arctan[cotan(i7nA/a:)]|7r- 1 . 



(53) 



The equal sign holds when \[/ has one of the values defined by (46) and the accuracy of the 
inequality (53) increases with decreasing a. Since, in general %w\p/a will not lie between 
and x the "reduced" property 



iTn///a=%inl//a±irl n , l n = 0, 1, 2, . . ., 
is introduced'and l n is such that \^\a is between and x. Thus 

| /(a^) | g x- 1 1 H(0) 1 1 arctan[cotan( W/«)]| , 



(54) 



and finally 



|/(a^)|^i|ff(0)||l-^/a|. 

4. Asymptotic Expansions 



(55) 



The derivation of an asymptotic expansion of the field U for large distances from the edge 
of the wedge is equivalent with the reduction of the integral (31) to an asymptotic series. Com- 
paratively simple is the case of a plane field incidence (fig. 5). Then for the exponential decay 
problem 

H(x)=exp(—yp cosh x). (56) 



</3 zTT-tp' f> = 7T 



*>=¥>' 




*P = 



Figure 5. 



ip z =7r+(p' 



[p=a 
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In this case F 2 (6) is represented as a Laplace transform with yp as parameter. Methods 
for the asymptotic development of integrals of this type are known [8, p. 45]. Since by (44) 
F 2 (S) is composed of two integrals of the form (45) it is sufficient to investigate 



I(a,\p) = — (2a)~ l sin (inl//a) [cosh (ttx/o)— cos (irip/a)] ] exp (— yp cosh x)dx, 

Jo 



(57) 



and yp has one of the values t± (</>±</> / ). The substitution 

cosh x=l+t (58) 



transforms (57) into 

/(«,$ = -exp (-yp) ( t K(t) exp (-ypt)dt, 

Jo 

with 



(59) 



K(t)=(2a)~ 1 sin (7r^/a)(2 + /)- 1/2 {cosh[T/a log (1 +/+ (/ 2 + 20"J- cosh (tt^/q:)}- 1 . (60) 

Thus (59) is represented as a Laplace transform. For the following, the discussion shall 
be restricted to case a, as described in section 3 and fig. 2. Cases b and c, as described in 
section 3, can be treated in the same manner. Now, K{t) as given by (60) represents a regular 
function in the vicinity of t = Q provided that <t>±(j)' t^7t (corresponding to yfi^Q by (46)). The 
singularities of K(t) lie on the negative axis and the smallest of their moduli is equal to 
Min [1+cos (0±</>')L 

The Taylor expansion of K{t) near the origin is 

K(t)=^A n (a^)t n } M<Min[l+cos(*±tf')], (61) 

71 = 

with 

Anfaf) = (n I)" 1 [d*K(t)l(dt*)] t . - (62) 

The first two coefficients are 

A (a^) = (2a)- 1 2- l/2 COinn(i7ryPla), 
A 1 (arf) = -(8a)- 1 2- 1/2 cotan(i7r^/a : ){l+2(7r/a) 2 [sm(§7r^/a)]- 2 }. (63) 

The asymptotic expansion of / (a x yp) for large yp is obtained by inserting the series (61) 
into (59) and integrating term by term [8, p. 48, theorem 6]. The result is 

I(a,f) exp(~7P)2^«(^)r(n+|)( 7 p)- w - 1/2 , -37r/2<arg 7 <3x/2. (64) 

n=() 

The coefficients A n tend to infinity when </>±<£' tends to ir (and hence \f/ to 0) and the expan- 
sion (64) becomes useless. This is the case when the configuration is such that one of the 
boundaries of geometric optics is approached. In this case the point / = becomes a singular 
point of K(t) and no expansion of the type (61) exists. In order to obtain an asymptotic ex- 
pansion which is valid over the whole region one proceeds as follows. The integrand of (57) 
at x = Q is equal to — (2a) _1 cotan (^inp/a). In the vicinity of ^=0 (the boundaries of geo- 
metric optics) 

-(2a)- 1 cotan(ix^/a) = -(7r^)- 1 + 0(iA) (65) 

and the singular part of (65) is independent of the angle a. Let the field U be now written 
in the form 

U=(U) a= ^ + [U~(U) a ^] (66) 

and by (28), (29), and (35) to (37) 

354 



r=ao- 2 .+^(</>-</OM^--^ (67) 

Now, (U) a -2* is given by (42), and by (44) 
Fi&±4>') — [F a (4>±<l>')] am u=I(a,T-4T*')^ 
The difference 

7(a,7r-*T*') -/(27r,7r-<^T0 / ) 



-/; 



(2a)" 1 sin f~ (x-0T0 r )] 



(4^-^08 (4*±4*0 



i/ zx fV / a-ta/\1 cosh(4x)-sin(4*±4*') 

COSll (TJC/a)— cos - (?r — * + *) 



exp (— 7p cosh ./') </.r (68) 



is by (65) such that the integrand of (68) vanishes identically when *±*' = 7r, and the 
asymptotic expansion of (68) is given as the difference of two series of the type (64) whose 
coefficients A n {a^— (/>+<£') — A n {2ir i ir— * + *') tend to zero when </>=F</>' tends to tt. Using (64) 
one arrives at the asymptotic expansion 

F t (*-*')-[^ 1 (*-0')]-lr±{^t(*+*')-[^(*+*')]--h}=S(7p) 

exp(-7p)Z)B re (7p)-"-« -3ir/2< arg7<37r/2. (69) 

The coefficients h u are defined by 

B ll =r(w+4){A w (a,T-*+*0-4n(2T,T— 

TM // (a,7T-</>-(/> , )-^»(27r,7r-(/>-(/) , )+^ /i (a,7r + (/> + *0 -A n (27r,7T+* + *0]} (70) 

With the A n given by (62), (63). It follows then from (67), (42), and (69) for the field excited 
by an incident plane field on a perfectly reflecting wedge of an angle a for case a (fig. 5) 

J7=exp [y P cos (*-*0]{l-4Erfc[(2 7 p) 1/2 cos (4*-4*')]} 

+ exp[ 7 pcos (0+0)]{l-4Erfc[(27p) 1/2 cos (4* + 4*')]}+S(*yp>, 
where S(yp) is asymptotically given by (69). It follows then for the wave problem 
£7=exp [ifcpcos (*-*0]{l-4Erfc[(2ifcp) 1 / 2 cos (4*-4*')l) 

Texp [ikp cos (0+*O]{l-4Erfc[(2ifcp) 1/2 cos (4*+4*0]}+#(^p), (71) 
where for large kp 

S ( ikp) - -exp (—ikp — iw/4:) X) B n i n (kp) - n ~^ (72) 

and these expressions are valid over the whole region of the wedge. The plus or minus sign 
has to be taken according as the acoustic or the electromagnetic case is considered. Obviously 
S represents the deviation of the wedge field from the half -plane field and the expansions (69) 
and (72) represent asymptotic expansions for this deviation term. If (yp) 1/2 cos (!</>+ §00 and 
(fcp) 1/2 cos(4^±4*') are large, i. e., for the far field at large distances from the boundaries of 
geometric optics one can replace the error function by its asymptotic expansion [9, p. 147], 



Erfc (z) -tt- 1 exp (-z 2 ) J2 (-l) n T{n+%)z- 2n ~\ -3tt/4< arg 2<3tt/4, 



n=Q 



(73) 



and the asymptotic expansion of U is of the same type as the one for S. On the boundary of 
the geometric reflection (<j>-\- 4,' = w ) , (71) becomes 
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£7=exp [ikp cos (<£-<£')]{l-i Erfc [(2ik P ) 1/2 cos (h<i>-W)]}Th exp (-ikp)+S(ik P ). (74) 
On the boundary of the geometric shadow ($— <l>' = w), 

[/=!exp (-0:p)=Fexp [(ifcp cos (0+</> / )]-{ 1 — i Erfc [(2ikp) 1/2 cos (|</>+i^)]}+S(^ P ). (75) 
Using (73) these results can be written as 

Z7=exp [ikp cos (</>— ^O] 1 ? i exp (-ikp) + 0[(kp)- 1/2 \, (74a) 

J7=Jexp (-ikp) + 0[(k P )-$]. (75a) 

These equations express the fact that the far field on the boundaries of geometric optics is 
independent of the angle of the wedge. 

5. Energy Radiated From a Hertz Dipole and From an Acoustic Point Source 
in the Presence of a Perfectly Reflecting Wedge or Corner 

Let the source Q(p', <t>', 0) be a Hertz dipole with its axis parallel to the edge. The func- 
tions/^) and H(x) are then 

f(0) = ( r 2 +2 2yi/2 e xp[-ik(r 2 + z 2 ) 1/2 ], (76) 

r=(p 2 + p ' 2 -2pp' cos0) 1/2 , 

H(x)=(p 2 +p' 2 +2pp' cosh x+z 2 )- 1 ' 2 exp[-^(p 2 +p ,2 +2pp / cosh x+z 2 ) 1/2 ]. (77) 

The function U given by (28) (taken with the minus sign) represents the vector potential 
A of the electromagnetic field. Its components are 

A x =A y =0, A Z =U. (78) 

With this relation the components E x , E Vi E z of the electric field vector E are known and 
the radiated energy- can be determined by Poyn tings method. The application of this method 
is very difficult in view of the complexity of the expression for U. It is therefore advisable 
to use the so called E. M. F. method for the determination of the radiation resistance R of 
the dipole [37, p. 457] and [33, p. 134]. R is given by 

R=-(L/J)Re(E 2 ) D . (79) 

Here L and J are the length and the current of the dipole respectively. He(E z ) D is the 
real part of the z component of the electric field at the location of the dipole Q(p', $', 0). It 
follows from (78) for the field, 

E t =-ic(^+¥Uy (80) 

<7is a constant (real ^>0) and is of no importance here since the "reduced" radiation resistance 
R/R will be determined. This is the quotient of the radiation resistance R for the complete 
system (dipole plus reflector) and the radiation resistance R of the dipole in free space (with- 
out the presence of the reflector). At first the contribution of F\ (which represents a set of 
dipoles consisting of the original source and a number of images) to E z is considered. From 
(30) and (76) 

U l =F l {<j>~4> f )-F l {<t>+^) 

V% g2 

= 2_) (?-|+2 2 )- 1/2 exp [(-• ^(rJ+3 2 ) 1/2 ]-2> s 2 +2 2 )~ I/2 exp [— ^(rf+s 2 ) 1 ' 2 ], (81) 

r 2 ,=p2-f-p'2_ 2 pp' cos (4>-4> / +2ap) 
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r 2 q =p i +p' i -2pp' cos (<j> + <t>' + 2aq) 

VTr+<t>-<t>'~\ r-x-<t>+4>'~\ 

Tir+^ + cp'"] r7r-0-^'"] 

2l= "L 2a J' ^ = |_~^~ 



It is then easily found that 



with 



d n = 2p' sin (an), d m =2p' sin (0' + am) 
™i = - [ Wa] i ^2= [ W a l , 



m 



=-[^i -[^ £ } 



j/(f) f ' sin f+r 1 ' cos f-r 3 sin f, 0(O)=2/3. 



(82) 



Re [ " ic (S i+Pf ■■'•) J = ~ 6 * 3 [,S ^^ ~ ,S - 9(Ww> ] (88) 



(84) 



The term m = in the sum (83) represents the contribution of the original dipole source. 
It follows then for the radiation resistance of the Hertz dipole in free space by (83) and (79) 



The contribution of F, to U is 



and it follows by (31 ) and (77) 



It = 2L/(3J)CP. 

r, = F,(<i>-4>')-F 2 (<t,+<t,') 



(85) 
(86) 



L V 2>^ 2 //>J Jo «Osh (Trxjot)— cos (7r 2 /a) 

,.„ v A sinp(7r-20ol sin p(7r+ 2^)1 

Jo 7 cosh (7r.//a)-cos ^ (x— 200 cosh (7rr/a)-cos - (tt+200 



,/./• 



r (87) 



It follows, then, for the quotient of the radiation resistance R of the dipole in the presence 
of the reflector and the radiation resistance R of the dipole in free space by eq (79), (83), (87), 
and (85). 



4a 



r n 2 1U2 "1 

«/ffo=3/2 2 ff(WJ- S <J{kd n ) - 

sinr-(7r-20')l 
-J o ,( 2 ^cosh0 



</ ( 2&p' cosh |j 

(irx/a)— COS (7r 2 /a) 



2 sin (tt 2 /c*) I , ,> "/ x </./; 

Jo cosh (ttx/o)— COS (7T 2 /< 



sm 



[>+ 2 ^] 



cosh (tjc/o)— cos - (x— 200 cosh (irx/a) — cos - (7T+20O 



<fa X- (88) 



The occurring parameters and #(f) are given in (84). Since </(f) vanishes for large f it is 
clear that if the distance of the dipole from the walls is large (p — >oo) only the term n— in the 
sum gives a contribution, and then 

lim JS/Z? = 3/20(0) =1. 
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The formula (88) is valid for an arbitrary a. But if a^>w, i. e., in the case of a wedge 
reflector (figs. 2 and 3), the two sums on the r. h. s. of (88) can be given explicitly. By (84) the 
limits n 1} n 2 , m X) m 2) of the sums in case a>7r are 



n ! = n 2 = rrii = m 2 = 
n l =n 2 = m l = 0, m 2 = 



for O<0'<K 
1 for £*<*'<£«. 



Thus 



MS 

z 



S »0WJ- S g(kd m )=2/3-g(2kp' sin </>'), if 0<*'<*t, 

= 2/3, if iT<*'<Ja. 

The analysis is exactly the same when the source Q{p ,<t>,0) is an acoustic point source. 
In this case the plus sign in (28) has to be taken. The function (/(f) in (84) which character- 
izes the coupling term between two Hertz dipoles (parallel to each other and the line con- 
necting their centers perpendicular to their axis) has to be replaced b} r the corresponding 
coupling term g x {f) for two acoustic point sources [30]. 



^i(f)=r 1 sinf. 



(89) 



In this case the corresponding formula for the relative radiation resistance as defined 
before is 



R/Ro= S <7i (fc/i.) + S (h (fcd m ) - (2a) - 
n=n\ m=m l 



2 sin (tt 2 , 



Jo 



(2k P ' cosh 
cosh (tx/cx)— cos (7r 2 a)' 



</i 



dx 



+ 



XX* 1 ' 00 * I) 



Sill 



[I fr- 2 *)] 



sin - 



[> +2 * r) ] 



h (inr/a) -cos 1"^ (tt-2p , )1 cosh (tt^)-cos j - (ir + 20')l 



cos 



dx 



(90) 

The definition of the n u n 2 , m u m 2 , d n , d m is the same as in (84), wliile (^(f) is given by 
(89). If c*y>T, then as before 

ni=n 2 =m 1 =m 2 =Q for <</> / <|tt, 

7ii =n 2 =m 1 =0, m 2 == — 1 for J'tt^^' < § a. 
and 



7*2 



hi , 

X 

m=m\ 



S ffi(fa/»)+ S flfi(WJ = l + (7i(2A:p / sin 00, 0O'<Jir, 



= 1, 



**'<>'<*«. 



6. Transient Solutions 

The results obtained in section 2 shall now be used to investigate the corresponding 
pulse problems. The majority of the contributions quoted at the beginning of this paper 
are based upon a direct attack on the time dependent wave equation. Integral transform 
methods were used in [38] (Kontorovich-Lebedev transform) and [40] (Fourier transform). 
In the cases considered here the application of Laplace's transform leads rapidly to the de- 
sired results. Let the excitation of the source be represented by a pulse starting at the time 
t=0. Let the pulse function be g(t), and g(t)=0 for £<0. Let the pulse function be such 
that it admits its representation by Laplace's integral formula: 
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;/(/) = (27ri)- 1 f C \i g(r) exp (-y T )dr 1 exp (yt)dy. 

Jc-ico Lt/0 J 



(91) 



In section 2, various expressions U for the exponential decay problem were given. This 
case corresponded according to (5) to the time dependency exp (yt), (put c=l). It can then 
be concluded thai for a pulse excitation g(t) as represented by (91) the expression for the field 
is represented by 



<P(t)=(2ri)- 1 i (j(r) exp (-7r) <It1 ( r exp (7/> 



/7- (92) 



The solution of the pulse problem is therefore given as the inverse Laplace transforma- 
tion of 

r* 00 

U g{r) exp (-7r>/r, (93) 

with respect to 7. 

6.1. Special Cases for g(t) 
For a unit stop function g(t) 1, the field (92) becomes 

$ (t) = (2tt0 - 1 1 Uy 1 exp (yf)dy. (94) 

%} c — i 00 

For a "Dime" pulse excitation imposed a1 the time £=£i^>0: 
g(t)=8(t—ti), one finds first 

g(r) exp (— 7r)dr= 6(7-/0 exp (— 7T)r/r = exp (—7/0. 
Jo Jo 

The field due to a "Dime" pulse excitation becomes then 

3> D {t-t l ) = (2iri)- 1 P ?7exp [y(t-U)]dy. 

J C-ia 

If (lie order of integration in (92) can be interchanged, it results in 
*((!)=(" g(T)I(2*i)- 1 \ U exp [y(t—T)]dy\dT. 



(95) 



(96) 



(97) 



Since by (96) the inner integral is the effect of a "Dime" pulse, excitation (07) can l>e 
written as 



*(*)=( 9(r)^ D (t-T)dr, 

Jo 



(98) 



i. e., the field due to an excitation g(f) of the source Q can be obtained by an integration proc- 
ess involving the field due to ;i "Dirac" pulse excitation. The general equations (96), (98) 
shall now be applied to the problem of the diffraction of arbitrary pulses (plane, cylindrical, 
and spherical) incident on a perfectly reflecting wedge or corner. The case of "Dirac" pulse 
excitation imposed at the time t=0 shall be considered first. The Laplace inversion of a 
function h(y) with respect to y shall be denoted by h(t), i. e., 

h(t) = (2iri) - 1 h exp (yt)dy. (99) 

%) c — i 00 

It follows from (96) and (28) for the field $ D due to a plane, cylindrical or spherical "Dirac" 

pulse: 

359 



* D (t)=F(<t>-<t>')^F(<p+<p') 

= J2f(<fi-<t>'+2an)T T, f(<l>+<t>'+2ctm)+\ H(i)[G(<f>-<j>')TG(<i>+<p')}(h; 

n = T)\ m—m\ Jo 



(100) 



where 



sin |~- (ir-0)"| sin |~- (tt+i9)1 

cosh (wx/a) — cos - (tt— 0) cosh (tx/o) — cos -(7r+0) 



<7(0) = -(2a)-N 



(101) 



The respective Laplace inversions /(0) and ff(#) are for a plane pulse: 

/(0) =6(/ + p cos 0) (102) 

H(.z)=d(t—p cosh x). (103) 

For a cylindrical pulse one has 

f(d)=2i<jr- l (t 2 -p 2 -p' 2 + 2pp' cos 0)- 1 ' 2 , for t>(p 2 + p' 2 -2pp' cos 0) 1 ' 2 , 

= 0, otherwise. (104) 

H(x)=2iir- 1 (f-p 2 -p f2 -2pp / cosh ./T 1 / 2 , for 2pp' cosh x<0 2 -p 2 -p' 2 , 

= 0, otherwise. (105) 

Finally, for a spherical pulse 

J(^) = ( r 2 + 32 ) -l/2 5 [^( r 2 + ^l/2] ( 1Q6 ) 

H(;x) = (p 2 + p' 2 +z 2 + 2pp'cosh x)- l <' 2 5[t-(p 2 + p' 2 + z 2 +2pp'cosh x) 1 ' 2 ]. (107) 

The Laplace inversions for /(0) and £TU) are obtained from (95) for the plane and for the 
spherical pulse, while [10, vol. 1, p. 277, formula 8] gives the inversion for the cylindrical 

pulse. It does not seem possible to express the integral I H(x)G(6)dx in a closed form in the 

case of a cylindrical pulse. To obtain the field due to a pulse excitation g(t), the equations 
(98) and (100) have to be used. For a plane and a spherical pulse the integration process in 
(98) can be performed immediately. The results are 

$(t) = Jl g[t + pcos(<t>-<t>'+2a>i)]3F S g[t+P cos (<p+<p'+2am)] 

n=ni m—mi 

/•cosh i=t/p 

+ (j(t-p cosh x)[G(4>-4>') TO^>+<t>')]dx 9 (108) 

for the plane pulse. 

m m-2 

+ (*' ( P 2 +p' 2 +s 2 +2pp' cosh x) - 1/2 g[t- (p 2 +p' 2 +s 2 +2pp' cosh x) " 2 ]{G(<p-<j>') T <?(<?+<?') }dx, 

(109) 
for the spherical pulse. 

The upper limit Xi of the integral is given by 

2pp'cosh Xi=t 2 — p 2 - P ' 2 -z 2 for (109), and cosh *i=*/p for (108). (110) 
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/;„ = [fr + p'*+r-2pp' cos ($-<!>' + 2 an)] 1 " 

(111) 
R n = [p 2 +p' 2 +z 2 -2pp' cos (0+<//+2<*m)] 1/2 . 

The other properties are given in (101). Again in the case of a wedge (c£>t), the prop- 
erties n iy n 2 , nil, "h are given in section 3. The representation of the cylindrical pulse field 
leads to a double integral. For the special case of a half plane, however, the integral occurring 
in (100) can be solved. In the following some special cases are considered in which the solution 
for the transient field become elementary. 

6.2. Special Cases 
a. Cylindrical Pulse and Half Plane 

The case of a line source Q parallel to the edge of a perfectly reflecting half plane is con- 
sidered. The regions I, II, III as described in section 3 are indicated in figure 6. For a half 
plane (a=2w), one obtains from (101), 

G(6)= -t~ 1 cos (id) cosh (&r)(cosh x+cosS)- 1 . (112) 

It follows then for a cylindrical "Dirac" pulse by (105), 



ip.ir-y 




<P 2 = 77><£' \,Q' (£,'-<?') 



FlGURK <>. 



f o "TI( X )G mx =-2i«-> cos (**)/; cosh (»*) ( '* *£££? * ) ~"' <fa - 

The upper limit x x of the integral is given by 2pp cosh x,i = t 2 — p 2 — p' 2 , and from (105) it 
is obvious that the integral is only different from zero when £> p J rp / . The substitution cosh 
x=\-\-v 2 , cosh {\x)dx—2 ll2 dv gives for the integral 



/; 



H{x)G(6)dx=-iir- 1 cos (|0)| sec (J0)|(* 2 -p 2 -p' 2 +2pp' cos0)- 1/2 , if t>p+p', 

= 0, otherwise. (113) 



It follows then from (100) by means of (35), (36), (37), (104), and (113) for the field due 
to a cylindrical " Dirac' ' pulse incident on a half plane 

<M*) = 22V-H(* 2 -Z? 2 )- 1/2 T(£ 2 -^^ for region I, 

= 2iT~ 1 { (£ 2 -R 2 )- 1/2 -M(* 2 -Ri)- 1/2 ± (P-Bl)- 112 }}, for region II, 

=%Tr- l [(t 2 -R\Y 112 ^ (t 2 -Rl)- l/2 ] 9 for region III, (114) 

where 

^=P 2 +p' 2 -2pp'cos(0-4/), i? 2 =p 2 +p' 2 -2 PP 'cos(<H-4>'). (H5) 

In (114), (t 2 —R 2 U2 )~ 1/2 =0 for t<Cp-\-p' when this expression is within the brackets; other- 
wise zero for £<i?!, 2 . The expression (114) is elementary in all the three regions. Former 
investigations [38] established $z>(0 as an infinite series involving Legendre functions. (It 
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can bo shown that this series can be summed and leads to the expressions in (114).) Fi'om 
(114) and (98) it follows for the field <1> due to an excitation g(t), 

<*>(*) =2wr- 1 / J ' h [(t-T)*-Rf\- l *g(T)dT* P " 2 [{t-TY-EV[- 1 'hj{T)dr 
-h£ "' [{{t-rY-RX)- l "T{{t-ry-I^-^]g{r)dr\, region I. 

$(0=2(T- 1 |£"" I [(<-r) 2 -i?r ,/2 .'7(r)'/r 

-hi' ' [((«-T)»-2??)- lfl ±(«-r)*-5S)- , «|ff(T)<iTJ»» region II. 

*)<j = *r- 1 J '[((t-r)'-^)" W: F(«-r) , -^)-wij(T)dT, region III. (116) 

b. Plane Unit Step Function Pulse and Wedge or Corner 

Now, g(t)=l for C>0, and g(t) = for /<0. Then the terms of the sums on the r. h. s. 
of (108) are either zero or one according as the argument of g is negative or positive. The 
well-known formula 

f [cosh (bx)-a]- 1 dx=2b- 1 (l-a 2 )- 1/2 arctan [(j^Y 2 tanh (Hoi 

gives for the integral in (108) by (101) 

Jo 0(6)d,= ^ a,ctan{i^ 2 -^L^cotan[£ (.-,)]} 

[ , + ( , 2 _ p 2 )1 / r / a+p ^ cotan|_-(x+^)j| 

for ty>p and equal to zero for t<Cp. Thus the field (108) for a unit step function excitation 
[19, 24] is elementary. 

c. Spherical Unit Step Function Pulse and Half Plane 

The center of the spherical pulse is Q (fig. 6). Again, g(f) = l. The excited held is 
given by (109) and G{6) is given by (112). The integral occuring in 109) is of the form 

P 1 (p 2 +p' 2 +2 P p cosh x)- 1/2 G(6)dx 

— -7T- 1 cos (|0) J cosh (ix) (cosh z+cos d)- l (p 2 +p' 2 +z 2 +2pp' cosh x)~ 1/2 dx 
= -2 1 / 2 tt- 1 cos (£0) ( l (1 + cos (9+^ 2 )- 1 [(p+p / ) 2 +2 2 +2ppV]- 1/2 ^, 

using the substitution cosh x=l+# 2 , cosh (^x)dx=2^ dv. The upper limit V\ is given by 
i) 1= =[t 2 — (p + p') 2 — z 2 ]K Again this integral is elementary. Each sum in (109) reduces to at 
most one elementary term (see the tables in section 3), and the resulting field [40] is elementary. 
It should be remembered that the previous investigations were based on the assumption 
that the velocity of the field is unity. The time parameter t has therefore to be replaced by 
ct when the velocity of the field is equal to c. The expression U for the time harmonic field as 
given in section 2 can easily be regained as the limiting case of a transient field $> due to an 
excitation g(t)=exp (iwt). 
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8. Appendix 
To prove 

irI v (a)K v (b) = I ' K [(a 2 +b 2 -2ab cos t) l/2 ] cos (vt) dt-sm (tv) 

Jo 

• f exp (-vt)K [(a 2 +b 2 +2ab cosh t) 1/2 ]dt, a<b, (Al) 

K v {a)K v (b)= [ K [(a 2 +b 2 +2ab cosh t) 1/2 ] cosh (vt) dt, (A2) 

K [(a 2 +b 2 -2ab cos <t>) l / 2 ] = 2w~ 1 [ K iv (a)K iv (b) cosh [v(w-\(j>\)}(h, 0^Re<^27r, (A3) 
one uses the formulas 

exp (z cos t) cos (vt) dt— sin (tv) exp (—z cosh t — vt) dt [9, p. 21], (A4) 
o ' Jo 

K v (z)=^[ x-'- 1 eiqp(-x~z*/x\dx, [9, p. 82], (A5) 

I v (a)K v Q>)=\ j" ^ 1 exp(-a;-^)/,Qa^)^, [9, p. 53]. (A6) 

8.1. Proof of (Al) 

Starting with (A6) and replacing the modified Bessel function in the integrand by (4), 
one obtains, interchanging the order of integration, 

rt \tt n\ l r / ,\ T P -i ( a 2 +b 2 -2ab cos t\ , 1 ,, 
TrI v (a)K v (b)=-\ cos (vt)\ x x exp( — x — — -\dx\dt 

1 • / n f" / ^rf OT i / a 2 ^b 2 +2abcosht\,l, 
--sm(7r^)J o ex P(-^)M Q x-'exp^-x — —^ -Jdxjdt. 

The inner integrals on the r. h.s. are known by (A5). This proves (Al). 

8.2. Proof of (A2) 
The relation 

K,(e) = iw sec(Tv)[I- p (z)-Uz)l (A7) 

together with (Al) gives (A2). 

8.3. Proof of (A3) 

Replace in (A2) v by iv and obtain by Fourier's inversion theorem 

K [(a 2 +b 2 +ab cosh t) 1/2 ]=2w~ 1 f " K tv (a)K iv (b) cos (vt) dv. (A8) 
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Replace here t by i (tt— </>) and obtain (A3). The restriction 0<Re <£<2tt is necessary to 
ensure convergence of the integral in (A3) at its upper limit, since for fixed real x and large 
positive v ^>x [9, p. 88, formula 19], 

K u (x) = (2ir) l i*(v 1 -&)- l,A e*v(^\*v\{smhv cosh-'(v/x)-(^-x')^+^+0(v~ l )\' (A9) 

8.4. Evaluation of S] Defined in (19) 



Using the formula 



* sin [(2N+l)y/2] 



one obtains instead of (19) 

Sl = (2tt) -> lim ( " A'„|7(p 2 +p' 2 -2 P p' cos x) ' ' 2 ] 

sin \j a (2JV+ 1) U+0)1 sin [^ (2A< + 1) (*-*)] 

Substitute in the first integral t(x + 6) = —2at and in the second jt(z— 6)=2aL It follows 
then that 



Ma 



«i 



=«r-«lim P 2 A' { 7 [pHp' 2 -2pp' cos (e+2atM)yn ^"f^^' dt, 
X-+0* Jt { sin t 



where ^i= -~ (tt+0) and t%=ir Or— 0)- Now *, is a Dirichlet integral and its value follows 
2a la 

immediately [2, p. 492] as represented in (22). 
Washington, D. C, January 14, 1958. 
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